Abstract. We present HST STIS observations of four 
quasar pairs with redshifts 0.84 < z cm < 1.56 and angular 
separation 2-3 arcmin corresponding to ~ 1— 1.5/ig Mpc 
transverse proper distance at z ~ 0.9. We study the dis- 
tribution of velocity differences between nearest neighbor 
Hi Lyman-a absorption lines detected in the spectra of 
adjacent QSOs in order to search for the possible corre- 
lation caused by the extent or the clustering properties 
of the structures traced by the absorption lines over such 
fS| ' a scale. The significance of the correlation signal is de- 
. termined by comparison with Monte-Carlo simulations of 
spectra with randomly distributed absorption lines. We 
find an excess of lines with a velocity separation smaller 
"> than = 500 km s" 1 significant at the 99.97% level. 
This clearly shows that the Lyman-a forest is correlated 
on scales larger than 1 Mpc at z ~ 1. However, out 
, of the 20 detected coincidences within this velocity bin, 
(NJ ; 12 have AV > 200 km s _1 . This probably reflects the 
fact that the scale probed by our observations is not re- 
r ~ l ' lated to the real size of individual absorbers but rather 
, to large scale correlation. Statistics are too small to con- 
ly-j ' elude about any difference between pairs separated by ei- 
£f) , ther 2 or 3 arcmin. A damped Lyman-a system is de- 
IT) tected at z abs = 1.2412 toward LBQS 0019-0145A with 
, log A^(Hi) ~ 20.5. From the absence of Znii absorption, 
' we derive a metallicity relative to solar [Zn/H] < —1.75. 
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1. Introduction 

Recent A^-body numerical simulations reproduce success- 
fully the global characteristics of the neutral hydrogen 
absorptions observed in quasar spectra, the so-called 



Lyman-a forest (Cen et al. 1994; 


Petitjean et al.| 1995, 


Hcrnc 


uist et al, 


1996; Zhang et al 


1995; Miicket et al. 


1996; 


Miralda-Escude et al. 


1996; Bond & Wadsley 1998). 



The absorptions arise from density inhomogeneitics in a 
smooth all-pervading intergalactic medium. Simulations 
show that the intergalatic gas traces the potential wells of 
the dark matter well at high redshift. It is therefore pos- 
sible to constrain the characteristics of the dark-matter 
density field from observation of the Lyman-a forest along 
a single line of sight (Croft et al. 2000). The addition of 
transverse information from observation of QSO pairs or 
more generally groups of quasars at small projected sep- 
aration in the sky will probably revolutionize this field 
in the next few years ( Petit jean| 1997). Indeed, inversion 



methods that have been recently implemented show that it 
is possible to recover the 3D topology of the dark-matter 
field using a dense network of lines of sight (Nusser & 
Haehnelt 1998, Pichon et al. 2001, Rollinde et al. 2001). 
After the early discovery of common absorptions in pairs 
of quasars (Shaver et al. 1982, Shaver & Robertson 1983, 
Weyman & Foltz 1983, Foltz et al. 1984), it was shown 
that the gaseous complexes giving rise to the absorp- 
tions should have large dimensions. In particular, stud- 
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ies of gravitationally lensed quasars (Smette et al. 1992, 
[Smette et al. 1995) yielded a lower limit of 100 h§Q kpc 
on the diameter of Lyman-a absorbers. Similar results 
were obtained from pairs of quasars with small separation 



separations have been 
(1998), |Dinshaw ct~aT 



(Bechtold et al. 1994; Dinshaw et al. 1994; Petitjean et al 
1998, p'Odorico et al.| 1998; pVIonier et al l 1998) 



Larger 



investigated 
(1998) 



by |Crotts fc Fang 



Monicr et alj (1999) and 



Williger et al. (2000). All studies conclude that absorp- 
tions are correlated on scales larger than 500 h$Q kpc. 
Unlike the case of QSO pairs with small angular separa- 
tions where the correlation can be explained by the fact 
that the lines of sight intercept the same absorber, the 
correlation for larger separations is certainly due to the 
clustering properties of distinct clouds. When observing 
triplets of quasars separated by 1 to 2 arcmin on the sky, 
corresponding to ^0.5 to 1 Mpc proper distance scales, 
both Crotts & Fang (1998) and Young et al. (2001) find 
statistically significant triple coincidences that they inter- 
pret as the presence of sheetlike structures along which 
inhomogeneous absorbers cluster. 

The number of such experiments is small however and it 
is important to increase the statistics. Here we present 
HST observations of four pairs of quasars. The ~ 2 and 
~ 3 arcmin angular separations between the two quasars 
of each pair probes scales between 1.0 and 1.5 ft,jj" Mpc 
proper distance at z ~ 1. This is where the transition be- 
tween individual halos and filamentary or sheet-like large 



scale structures is expected (Miicket et al. 1996, Charlton 
1997). 

We describe the observations in Sect. 2 and comment on 
individual metal line systems in Sect. 3. Correlations be- 
tween metal line and Lyman-a systems are respectively 
discussed in Sect. 4 and 5. Conclusions are drawn in 
Sect. 6. 
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Fig. 1. HST STIS spectra of the eight observed quasars. The tick marks indicate significant (>4cr) absorption features. 
The vertical dashed lines mark the red limit of the wavelength range for the selection of the Lyman-a lines. This limit 
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2. Observations 



Table 1. List of observed QSOs 



Object name 



A8 



Az 



< S > 



LBQS00f9-0145A 1.59 

LBQS0019-0145B 1.04 3.3 0.70 - 1.01 1640 
Q0035-3518 1.20 

Q0035-3520 1.52 3.4 0.69 - 1.17 1710 

Q0037-3545 1.10 

Q0037-3544 0.84 1.7 0.59 - 0.82 810 

PC1320+4755A 1.56 

PC1320+4755B 1.11 1.9 0.62 - 1.08 940 
a Angular separation on the sky in arcmin 
6 Redshift range over which coincidences were searched for 
c Mean proper distance in kpc between lines of sight in the 
redshift range (q = 0.5, A = 0, H = 50 km s -1 Mpc -1 ) 

Observations were carried out on the Hubble Space Tele- 
scope using the Space Telescope Imaging Spectrograph 
(STIS) with the G230L grating and the Near-UV-MAMA 
detector. This configuration yields a mean spectral reso- 
lution of R = 700 (FWHM = 3.4 A at A = 2374 A) and a 
wavelength coverage from 1570 A to 3180 A. The obser- 
vations were reduced at the Goddard Space Flight Center 
with the STIS Investigation Definition Team (IDT) ver- 
sion of CALSTIS ( |Lindler| 1998). Standard reduction and 
calibration were used. Special care was taken to determine 
accurately the background due to the sky and the dark 
current. The zero point of the wavelength scale for indi- 
vidual exposures was determined requiring the Galatic in- 
terstellar absorptions to occur at rest. The correction can 
always be performed because the Galactic Mgn doublet 
is well detected in every single spectrum. When the Fen 
lines were also detected we checked that the dispersion in 
the zero point is smaller than the spectral resolution. The 
resulting spectra are shown in Fig. 1. The quasar contin- 
uum was fitted with Gaussian profiles for emission lines 
and simple cubic splines in regions between emission lines. 
The best fit was found by varying the position of the con- 
trol points of the cubic splines. The resulting continuum 
was slightly manually adjusted in regions that were poorly 
fitted as for example near broad emission lines and Lyman 
limits. The detection of absorption lines in the normalized 
spectrum was performed by filtering each spectrum to im- 
prove the contrast between lines and noise. This has been 
performed by successively applying a wavelet filter and an 
"upgraded" median filter to the spectrum. 
We used for the wavelet filter B3-spline scaling functions. 
This filter selects the wavelength scales corresponding to 
the width of the absorption lines (see panel b of Fig. ^). 
Pixels from the wavelet filtered spectra are sorted in in- 
creasing order keeping trace of the pixel permutations. 
The distribution of the pixel values is shown in panel (c) 
of Fig. A lower limit of the level of noise in each pixel 
can be estimated (dashed line) and substracted to the real 



pixel value (panel d in Fig. 2). Pixels are then reordered 
and the resulting spectrum is shown in panel (e) of Fig. |[ 
This filtered spectrum is used to define regions where pos- 
sible absorption lines are present using a threshold defined 
so that no line is lost in the next step. We then compute 
in the original spectrum the equivalent width and the as- 
sociated noise over each of these regions and select only 
those with an equivalent width to noise ratio larger than 

2. The corresponding regions are then fitted with a Voigt 
profile fitting program to derive the position of absorption 
features. This software makes a x 2 minimisation in each 
region adding lines until the reduced % 2 reaches a value 
lower than or equal to 1. 

For each absorption feature, we calculate the equivalent 
width in windows centered on the minimum of the line 
and of widths an increasing number of pixels. The signal- 
to-noise ratio computed from the noise spectrum is plotted 
as a function of the distance to the central pixel (panel f of 
Fig. 2). The S/N ratio of the line was taken at the maxi- 
mum of the curve and the equivalent width was computed 
by integrating the fitted profile. The lines with S/N ratio 
greater than 4 are listed with their identification in Tables 
2 to 5. In these tables, uncertain positions or identifica- 
tions are indicated by a colon. 

3. Comments on individual metal line systems 

In this Section we comment on the intervening metal line 
systems identified in the spectra. To refer to an absorption 
feature, we use the numbering given in Tables 2 to 5. 



3.1. LBQS 0019-0145A 



1.59 



The HST data on this quasar have been comple- 
mented with a high resolution (R ~ 40000) UVES spec- 
trum covering the wavelength ranges 3900-5200 A and 
5450—9300 A. The exposure time was two hours. 

3.1.1. z abs = 0.6514 

The strong absorption feature #13 at A2559.23 cannot be 
explained by Lyman-/? at z a b s = 1.4976 alone because the 
corresponding equivalent width is too large. As there is 
a strong Civ system at z a b s = 0.65142 in the spectrum 
of LBQS 0019-0145B, we tentatively note that the addi- 
tional absorption could be due to Civ at z a b s = 0.6514. 
The corresponding Lyman-a absorption is unfortunately 
redshifted below the Lyman limit of the z a b s = 1.4976 
system. MgnA2803 at this redshift is not detected in the 
UVES spectrum down to w bs = 35 mA. 

3.1.2. z abs = 0.6953 

MgnA2796 (w r ~0.8 A), FeiiA2382 (to r ~0.5 A), MgiA2852 
and Call absorptions are detected at this redshift in 
the UVES spectrum in addition to AlnA1670 (line #28) 
which is detected in the HST spectrum. Unfortunately, 
the HiA1215 line is redshifted at A = 2060 A, below the 
Lyman limit of the z a hs = f -4976 system. 
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Table 2. Line list for the pair LBQS0019-0145AB 



Table 3. Line list for the pair LBQS0035-3518 & Q0035- 
3520 
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30 


2900 


20 





89 


Lya 


1 


38568 
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The profile of the Mgn and Fen absorptions consists of two 
main absorption features approximately 100 km s" 1 apart 
and separated by a sharp drop in optical depth near the 
center (see Fig. ||). As described by Bond et al. (2001), 
and despite the moderate MgnA2796 equivalent width, 
this profile may indicate that the line of sight intercepts a 
superwind. 
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3.1.3. z abs = 1.2412 

Strong HiA1215 (line #21), HiA1025, Sin, SimA1206, Siiv, 
Nv and CnA1334 absorptions associated with this sys- 
tem are detected in the HST spectrum. Additional Mgn, 
MgiA2852, Fen and A1iiiA1854, 1862 absorptions are seen 
in the UVES spectrum (see Fig. ||). Despite the low resolu- 
tion of the HST spectrum, damped wings are clearly seen 
in the case of the Lyman-a line. The simultaneous fit of the 
Lyman-a, and Lyman-/3 absorptions gives log N (Hi)^20.5 
(see Fig. 4). The observed Fen and Mgn lines are heavily 
saturated so that only a lower limit on the column den- 
sities can be derived, log N(Fen & Mgn)>15. Znll is not 
detected with log A r (Znll) < 11.4 implying that metallic- 
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Table 4. Line list for the pair Q0037-3544 & Q0037-3545 Table 5. Line list for the pair PC1320+4755AB 
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Fig. 3. Metal absorption lines at z — 0.6953 (left column) 
and z — 1.2412 (right column, damped Lyman-a system) 
in the LBQS 0019-0145A UVES spectrum. 



ity relative to solar, [Zn/H] < —1.75, is one of the smallest 
observed at low redshift (see e.g. Lcdoux ct al. 2001). Note 



that the large column density found for Mgn implies that 
the neutral part of the system cannot account for all the 
Mgn we see. 
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As for the system at z a b s = 0.6953, the profiles of the ab- 
sorptions consist of two main components separated by 
about 50 km s . This again may indicate that the line of 
sight intercepts a superwind. 

3.1.4. z abs = 1-3088 

The system is unambiguously identified by sharp 
MgnA2796, 2803 absorptions (w r ~0.07 A) detected in the 
UVES spectrum. The Lyman-a absorption is blended with 
MgnA2796 from the interstellar medium. 

3.1.5. z abs = 1-4976 

This strong system is detected by HiA1215, 1025, 972, 949, 
937 and CiiiA977 absorptions. It is at the origin of the 
Lyman limit at A — 2270 A and therefore log N(Ri) > 18. 
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Fig. 4. Lyman-a and Lyman-/3 absorptions centered at 
z abs = 1.2412 toward LBQS 0019-0145A. The continuous 
curve is the best fit to the two Hi absorptions obtained 
with log AT(Hl) ~ 20.5. 

A possible MgnA2796 absorption line is present in the 
UVES spectrum with w bs ~ 25 mA. 

3.2. LBQS 0019-0145B z cm = 1.04 

Only one metal line system along this line of sight is de- 
tected at z abs = 0.6513 by HiA1215, CivA1548, 1550 and 
SiivA1393, 1402 absorptions. 

3.3. Q 0035-3518 z cm =1.20 

3.3.1. z abs = 1-088 

Strong HiA972, 1025, 1215 and SinA1260 absorptions are 
detected in this system. The latter line is probably blended 
with an other Lyman-a line however. Moreover, an absorp- 
tion feature detected at the 2.5cr level is observed at the 
expected position of CnA1334 at A~2785 A. 

3.3.2. z abs = 1.1954 

This system is at slightly larger redshift than the quasar 
(~ +1200 km s _1 ). It shows strong associated Ovi absorp- 
tion. NvA1238 is detected at the 2.5er level. 

3.4. Q 0035-3520 z om =1.52 

The four metal line systems detected in this quasar are 
all within 3500 km s _1 from the emission redshift of the 
quasar. As they are approximately at the same redshift, 
their Lyman-/3 and Ovi absorptions are blended. The re- 
sulting blend corresponds to the four strong features seen 
at A-2600 A in Fig. |. 



3.4.1. zabs = 1-4927 

The Lyman-/3 line of this system is observed at A~2550 A, 
but is below the 4<r detection limit. As no other metal line 
is detected, the CiiiA977 identification is tentative and 
the feature at A2437.38 could be a blend of Lyman-7 at 
z a bs = 1.50621 and Lyman-a at z a b s = 1.00492. 

3.4.2. Zabs = 1.50362 

This system consists of HiA949, 972, 1025, 1215, Ovi and 
CiiiA977 absorptions. A feature is also observed at the ex- 
pected wavelength of the associated SimA1206 absorption 
(A-3020 A). 

3.4.3. zabs = 1-5118 

HiA1025, 1215, Ovi, CiiiA977 and SimA1206 are observed 
at this redshift. 

3.4.4. Zabs = 1-5158 

This system has all the characteristics of an associated 
system with very strong Nv and Ovi lines. Therefore the 
Nill and Cm identifications are tentative. 

3.5. Q 0037-3544 zcm = 0.84 

3.5.1. zabs = 0.5922 

Strong HiA1215 and Civ absorptions are detected. 
SiivA1393, 1402 are redshifted at A2219.4 and 
2231.31 A (lines #8 and #9). Given the strengths 
of the corresponding absorption features, it is most prob- 
able that the two Siiv lines are blended with additional 
Lyman-a lines. An absorption feature is detected at the 
expected position of SiinA1206 (A~1920 A) but is below 
the 4cr threshold. 

3.5.2. zabs = 0.6961 

This system consists of a Civ doublet detected outside the 
Lyman-a forest. The Lyman-a absorption is observed at 
wavelength 2062 A, but is below the 4er threshold. 

3.5.3. zabs = 0.8364 

HiA1215, SiinA1206 and Civ absorptions are observed at 
this redshift. Moreover, a feature is observed at the ex- 
pected position of SinA1260 (A-2316 A). 

3.6. Q 0037-3545 z cm =1.10 

We do not identify any metal line system along this line 
of sight. 

3.7. PC 1320+4755 A z cm =1.56 
3.7.1. zabs = 1-0762 

Strong HiA949, 972, 1025, 1215, CmA977, Siiv, SinA1260 
and CnA1334 absorptions are detected at this red- 
shift. Additional features are observed at the expected 
wavelengths of the associated Niii989 (A~2055 A) and 
SinA1190, 1993 (A-2470 & 2480 A). This system is at the 
origin of the Lyman limit seen at ~ 1890 A and therefore 
log AT(Hl) > 18. 
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Table 6. List of metal lines for the four pairs Table 7. List of Lyman-a lines for the four pairs 
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3.7.2. z abs = 1.4329 

Strong Hi (from Lyman-a to Lyman<) , Ovi, SimA1206 
and SinA1260 absorptions are observed. From the partial 
Lyman limit seen at ~ 2210 A, we derive log AT(Hi) ~ 17.3. 

3.8. PC 1320+4755B z cm =1.11 

Only one metal line system is detected along this line of 
sight at z abs = 0.9270. It consists of HlA1215, SimA1206, 
Siiv and Civ absorptions. A feature is observed at the 
position of the associated SinA1260 (A— 2430 A). 

4. Correlation of metal line systems 

We summarize in Table 6 the metal line systems seen along 
one line of sight for which the corresponding absorption 
along the adjacent line of sight could be observed. In case 
of blending, w bs,min and w bs,max are the minimum and 
maximum equivalent width of the transition. These values 
are computed using consistency arguments relating the 
equivalent widths of lines observed in the same system 
(see examples in Sect. 5.1). The 4cr detection limits are 
also indicated. 

There are only 6 metal line systems which have all 
w r (Lyman-a) > 1 A. Apart from the Civ system at 
z abs = 0.65 toward LBQS 0019-0145B which may 
have a coincident absorption in LBQS 0019— 0145A (see 
Sect. 3.1.1), none of the other systems are detected 
along the adjacent line of sight down to a 4a limit of 
w r (Lyman-a) < 0.40 A. As the system at z a b s = 1.19 to- 
ward Q 0035—3518 could be associated with the quasar, 
this means that out of 5 intervening metal line systems 
with w r (Lyman-a) > 1 A, only one is present in the two 
lines of sight. 

Correlation of Civ systems however has been claimed on 
large scales at high redshift (e.g. Williger et al. 1996). 
Moreover, of the five w r > 0.4 A Lyman-a systems seen 
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calculated using the width of the line detected along 



cither line of sight. 

at the same z ~ 2 redshift in the three spectra of KP 76, 77 
and 78 (triple hits over 2-3 arcmin separations), two show 
associated Civ although Civ is seen only in about one 
w T > 0.4 A Lyman-a system out of ten (Crotts & Fang 
1998). All this may indicate that the transverse clustering 
of Civ systems is less pronounced at z ~ 1 than at z ~ 2 
(see also D'Odorico et al. 2002). 

5. Correlation in the Lyman-a forest 

5.1. The Lyman-a line list 

From the line lists obtained as described in Sect. 2 and 
given in Tables 2 to 5, we have extracted for each pair 
of QSOs, a master line-list of Lyman-a lines based on 
several criteria: (i) we include all isolated lines when no 
other identification is found; (ii) the lines must be at more 
than 3000 km s _1 blueward of the two Lyman-a emission 
lines; (iii) we use physical consistency arguments to infer 
the presence of Lyman-a lines blended with metal lines 
(for example CivA1548 cannot be weaker than CivA1550); 
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only limits on the equivalent width can be inferred this 
way; (iv) we impose some equivalent width threshold. The 
Lyman-a line list is summarized in Table 7. The columns 
correspond to: #1 and #6 line number in the spectrum 
(see Fig. 0); #2 and #7 Lyman-a redshift; #3:#4 and 
#8:#9 the range in equivalent widths in case of blending 
(see below); #5 and #10 the 4cr equivalent width detection 
limit using the width of the line detected along either line 
of sight. 

In LBQS 0019-0145A, absorption #3 cannot only 
be OviA1037 as it has w b s — 2.10 A whereas 
^obs(OVI1031) = 1.64 A. The hidden Lyman-a line has 
0.28 < w bs < 0.49 A corresponding to the Ovi doublet 
ratio ranging from 1 to 2. 

Absorption #8 coincides in redshift with Cm at 
^abs = 1.49672. However, we consider this identification 
unlikely. Indeed, the line is quite strong (w b s = 1-12 A) 
even though the UVES and STIS spectra do not show any 
other metal lines at this redshift except for a very weak 
MgnA2796 line system and a 2.5a feature shifted by 1.5 A 
from the expected position of SimA1206. Therefore, we 
identify this line as Lyman-a at z a b s = 1.00660. 
In LBQS 0019-0145B, the limits on the equivalent 
width of the Lyman-a line blended with SiivA1402 at 
^abs = 0.65381 in feature #12 are derived applying the 
doublet ratio to the SiivA1393 equivalent width. 
In Q 0035—3518, we estimate the equivalent width of 
Lyman-7 (line #2) from the associated Lyman-a and 
Lyman- (3 absorptions. 

In Q 0035-3520, for feature #7, we assume that CmA977 
at z a b s = 1.49469 contributes very little since the associ- 
ated Lyman-a absorption is weak and no other metal line 
is detected; limits on the equivalent width of the Lyman- 
a line at z a b s = 1.00497 come from limits on Lyman-7 
at z a bs = 1.50621 derived from the Lyman-a absorption. 
Similarly, we have constrained the equivalent width of 
Lyman-/? at z a b s = 1.47617 in the absorption feature #11 
from the associated Lyman-a line. 

In PC 1320+4755A, the limits on the equivalent width 
of Lyman-<5 (feature #1) and Lyman-7 (feature #2) at 
z a bs — 1.07575 come from their associated Lyman- (3 and 
Lyman-a absorptions. For feature #3 we assumed for 
the equivalent width of CinA977 a conservative limit of 
1.5 A from the associated CnA1334 line. Finally, lower 
and upper limits on the equivalent width of OviA1031 at 
z abs = 1.07726 (feature #5) and FeiiA2344 at z abs = 
(feature #11) come from OviA1037 which has an equiva- 
lent width w o b s ^0.8 A, below the 4tr detection limit, and 
FeiiA2600 respectively. 

Note that the number of metal lines which could be 
misidentified as Lyman-a lines is expected to be small. 
Indeed, they can be neither Mg II nor Fe 11 lines because 
the strongest potential lines of these species are redshifted 
beyond the wavelength of the Lyman-a line with the high- 
est redshift in our sample. They cannot be lines too close 
in wavelength to Lyman-a as our careful procedure would 



have identified them. The only possibility is that some 
isolated Al nA1670 or C ivA1548 (with C ivA1550 not de- 
tectable) lines could be present in the wavelength range 
where we search for coincidences. The number of C IV and 
Mg II systems with w T > 0.3 A at z - 0.3 is 0.87 and 0.75 
per unit redshift respectively (Bergeron et al. 1994, Boisse 
et al. 1992). We consider that half of the Mg 11 systems 
have an associated Al 11 line and half of the C IV systems 
would have only C IVA1548 detected. Moreover, only less 
than two-third of these systems have w r > 0.4 A which is 
the usual 4 a limit of our spectra. As our survey samples 
a redshift interval of Az = 2.96 (considering six lines of 
sights, see Table 1), the expected number of misidentified 
lines is of the order of 1 to 2. This is to be compared to the 
46 Lyman-a lines with w r > 0.4 A we detect. Note that 
the latter number is consistent with the number of lines 
detected in the HST Key-program (Jannuzi et al. 1998). 

5.2. Correlation 

From the master line list, we selected lines with w r >0.3 A 
and applied the N earest-Neighbor method as described in 
e.g. Young et alj (2001) to estimate the level of correla- 
tion between absorptions detected along adjacent lines of 
sight. In this method, there is no a priori velocity separa- 
tion limit in the definition of a coincidence. A couple of 
lines along two different lines of sight is declared to be a 
coincidence if each of the lines is the nearest neighbor of 
the other. Note that the procedure underestimate the clus- 
tering signal as it does not take into account the difference 
of S/N ratio along the two lines of sight. 
The distributions of velocity separations (|Ai>|) between 
the two lines in a coincidence are plotted in Fig. || for 
the pairs with angular separation ~2 arcmin (Panel a), 
3 arcmin (Panel b) and the complete sample (Panel c). To 
im prove statist ics, we have added to our results, the data 
of |Young et al. (2001) on two additional pairs separated 
by ~2 arcmin and ^3 arcmin. The results are plotted in 
Fig. H for the pairs with an angular separation ~2 (Panel 
d), 3 arcmin (Panel e) and for all the pairs (Panel f). 
To estimate the excess of correlation with respect to ran- 
domly placed absorption lines, we produced 100000 sim- 
ulated master line lists drawn from a population of ran- 
domly redshifted lines, taking the same number of lines 
and the same wavelength range as in the observed spectra. 
Results of applying the same method to the simulated line 
lists are given as dotted lines in Fig. ||. The error bars in 
the figure correspond to the rms of the values found in the 
simulation. As the corresponding distribution is not Gaus- 
sian, we indicate in the following the probability that the 
observed number of coincidences occurs in the simulated 
population. 

We detect 2 and 4 coincidences with \Av\ smaller than 
500 km s _1 in the two pairs separated by ^2 arcmin and 
^3 arcmin respectively. The fact that the number of coin- 
cidences is smaller in the closest pairs is not statistically 



10 



Bastien Aracil et al.: HST STIS observations of four QSO pairs 



10 



5- 



10 




20 



5- 



10- 



1000 



2000 



3000 



(a) 



(b) 



4000 





( C ) £ 














•■"4- J H 4n I r^ V .\ 



(d) - 



1000 



2000 



3000 



4000 



(e) - 



(f) - 



in 



+ 



1000 2000 3000 

Relative velocity (km s" 1 ) 



4000 



Fig. 5. Number of coincidences versus the velocity separa- 
tion | Aw | between the two lines for rest equivalent width 
threshold of w r > 0.3 A. Panel (a): for pairs with sep- 
aration - 2 arcmin, Q 0037-3544 & Q 0037-3545 and 
PC 1320+4755A&B; Panel (b): for pairs with separation 
- 3 arcmin, LBQS 0019-0145A&B and Q 0035-3518 & 
Q 0035—3520; Panel (c): for the four pairs. Panels (d),(e) 
and (f) are the same as, respectively, Panels (a), (b) and 
(c) but after adding data of Young et al. 2001. Dotted 
lines correspond to the expected number of coincidences 
from a randomly placed population of lines. 
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Fig. 6. Number of coincidences versus the velocity sepa- 
ration | Aw| between the two lines for the complete sample. 
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7. Velocity distribution of the 20 coincidences with 
< 500 km s _1 . The dotted line indicates the expected 



number of coincidences from the Monte-Carlo simulations. 



significant. When the complete sample is used (Panel c of 
Fig.[|), the excess in the first bin (|Aw|<500 km s _1 ) is sig- 
nificant at the 99.20% level. To increase the statistics, we 



have added to our sample data from Young et al. (2001). 
The total number of coincidences with |Aw|<500 km s _1 
and w r > 0.3 A is increased to 6, 8 and 14 for the pairs sep- 
arated by ~2 arcmin, ^3 arcmin and the complete sample. 
The corresponding excesses relative to simulations of ran- 
domly placed lines are significant at the 95.57%, 99.92% 
and 99.97% levels respectively. The excess is about the 
same when no threshold is applied to the equivalent width 
(see Fig. [|). In that case, the number of coincidences in 



the first bin is 20 and the excess is detected at the 4<r level. 

This clearly shows that the Lyman-a forest is correlated 
on scales larger than 1 h$Q Mpc proper at z ~ 1. However, 
we should note that, in the complete sample, 12 of the 20 
coincidences with Aw| smaller than 500 km s _1 actually 
have | Aw | greater than 200 km s _1 . This is also the case 
for 8 of the 14 coincidences in the wj r > 0.3 A sample. 
These velocity differences are probably related to peculiar 
velocities of different objects and reveal that the scale we 
probe is not related to the real size of individual absorbers. 
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Fig. 8. Equivalent width observed along one line of sight 
versus the one observed along the second line of sight for 
the 20 coincidences with velocity difference smaller than 
500 km s _1 . 

In Fig. H we plot the equivalent width observed along one 
line of sight versus the equivalent width observed along 
the second line of sight for the 20 pairs with velocity dif- 
ferences smaller than 500 km s _1 . There is no clear trend 
in the plot. This is not really surprizing as we do not ex- 
pect the absorption strengths to be correlated at such a 
separation neither in the case of a common absorber nor 
in the case of independent objects. 

6. Conclusion 

We have searched for coincidences of Lyman-a absorbers 
at z < 1 along the lines of sight toward four pairs of 
quasars with separations 2—3 arcmin, observed with HST 
STIS. Using the Nearest-Neighbor statistics, we have con- 
structed the distribution of the velocity difference between 
absorption lines detected along two adjacent lines of sight. 
We have compared this observed distribution to that de- 
rived from Monte-Carlo simulations placing at random the 
same number of lines in the same wavelength ranges as in 
the observations. For lines with w T > 0.3 A , we find an ex- 
cess of coincidences with velocity separations smaller than 
500 km s _1 significant at the 99.2% level. Combining our 
data with those in the literature (Young et al. 2001), the 
excess relative to simulations is significant at the 99.97% 
level for lines with w T > 0.3 A and at the 99.98% level if 
no condition on the equivalent width is imposed. 
The result is consistent with the findings of similar stud- 



This is not a systematic effect related to the low spec- 
tral resolution of our data. It could be related to the in- 
crease of peculiar velocities with decreasing redshift for 
comparable spatial scales. This conclusion is strength- 
ened by the fact that the transverse correlation, x K 
£ 4 (1 - Fi(Ai)) x (1 - F 2 (X l )) where F^K) and F 2 (A 4 ) 
are the two QSO normalized fluxes at wavelength A^, at 
z ~ 1 measured on our spectra is very small, \x\ < 0.01, 
whereas it is of the order of 0.2 at z ~ 2 for the same sep- 
aration (Rollinde et al. in preparation, see also McDonald 
2000, 2001 and Viel et al. 2001). This can be explained 
by the fact that large values of \ are due to coincidences 
with velocity separations smaller than the spectral reso- 
lution. In conclusion, evolution from high to low redshift 
is not seen in the level of correlation but rather in the 
velocity difference between lines of sight which increases 
with decreasing redshift. Simulations have shown that ab- 
sorption lines with a given column density correspond to 
higher overdensities at low-redshift compared to high red- 
shift and most of the Lyman-a forest at low redshift is 



to be revealed by very weak lines (Ricdiger et al. 1998; 
Thcuns et al] 1998; Penton et al. 2000). Therefore, de- 
tailed study of the cosmological evolution of the Lyman-a 
forest will only be possible when the sensitivity of the in- 
struments will be high enough to routinely detect lines 
with w T < 0.01 A in the UV. 
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Fig. 2. Illustration of the procedure used to detect absorption lines. Panel (a) represents the original normalized 
spectrum which is wavelet filtered to select scales characteristic of the absorption lines. The filtered spectrum is shown 
in panel b. All pixels in the spectrum are sorted according to their value (panel c). Low values correspond to absorption 
and high values to spikes; the intermediate values are dominated by noise which is fitted by the dashed line in panel 
(c). This level of noise is subtracted to each pixel to give the new distribution of panel (d) and pixels are reordered. 
The resulting spectrum is shown in panel (e). This spectrum is used only to select the regions where absorption lines 
are detected. The equivalent width of each absorption feature is measured in the original spectrum. For each line, the 
signal to noise ratio is plotted as a function of the distance to the central pixel (panel f). The S/N ratio of the line is 
taken at the maximum of the curve. 



